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ABSTRACT 

We present Roche tomograms of the G5-G8 IV/ V secondary star in the long- 
period cataclysmic variable BV Cen reconstructed from MIKE echelle data taken on 
the Magellan Clay 6.5-m telescope. The tomograms show the presence of a number 
of large, cool starspots on BV Cen for the first time. In particular, we find a large 
high-latitude spot which is deflected from the rotational axis in the same direction 
as seen on the K3-K5 IV/V secondary star in the cataclysmic variable AE Aqr. BV 
Cen also shows a similar relative paucity of spots at latitudes between 40-50° when 
compared with AE Aqr. Furthermore, we find evidence for an increased spot coverage 
around longitudes facing the white dwarf which supports models invoking starspots 
at the L\ point to explain the low-states observed in some cataclysmic variables. In 
total, we estimate that some 25 per cent of the northern hemisphere of BV Cen is 
spotted. 

We also find evidence for a faint, narrow, transient emission line with charac- 
teristics reminiscent of the peculiar low-velocity emission features observed in some 
outbursting dwarf novae. We interpret this feature as a slingshot prominence from 
the secondary star and derive a maximum source size of 75,000 km and a minimum 
altitude of 160,000 km above the orbital plane for the prominence. 

The entropy landscape technique was applied to determine the system parameters 
of BV Cen. We find Mi = 1.18 ±%H M , M 2 = 1.05 ±g;g M Q and an orbital 
inclination of i = 53°± 4° at an optimal systemic velocity of 7 = -22.3 km s _1 . Finally, 
we also report on the previously unknown binarity of the G5IV star HD 220492. 

Key words: stars: novae, cataclysmic variables - stars: spots - stars: late-type - 
stars: imaging - stars: individual: BV Cen - techniques: spectroscopic 



1 INTRODUCTION 

Cataclysmic variables (CVs) are short period binary sys- 
tems in which a (typically) late main-sequence star (the 
secondary) transfers material via Roche-lobe overflow to a 
whi te dwarf prima ry star. For an excellent review of CVs, 
see IWarnerl (|l995h . Although CVs are largely observed to 
study the fundamental astrophysical process of accretion, 
it is the Roche-lobe filling secondary stars themselves that 
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are key to our understanding of the origin, evolution and 
behaviour of this class of interacting binary. 

In particular, the magnetic field of the secondary star 
is thought to play a crucial role in the evolution of CVs 
driving C Vs to shorter orbital per i ods through magnetic 



braking (e.g.lKraftlll967l,lMestellll968l . lSpruit fc Ritte: 



lagnenc 

iRappaport. Verbunt fc JossHl983T ). Furthermore, the transi- 
tion of the secondary star to a fully convective state and the 
supposed shutdown of magnetic activity that this transition 
brings is also invoked to explain the period gap - the dearth 
of CVs with orbital periods between ~ 2-3 hours. On more 
immediate timescales, magnetic activity on the secondary 
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stars are also thought to explain variations in CV orbital 
periods, mean bright nesses, mean outburst durations and 
outburst shape s (e.g. lBianchinlll990l : iRichman et al I I 1994 
lAk et alJlioOll ). 

That the secondary star and its magnetic field should 
have such a large and wide-ranging impact on CV prop- 
erties should come as no surprise - the secondary star es- 
sentially acts as the fuel reserve that powers these binaries. 
Therefore, an understanding of the magnetic field proper- 
ties of the secondary stars in CVs (e.g. spot sizes, distri- 
butions and their variation with time) is crucial if we are 
to understand the behaviour of these binaries. In addition, 
detailed studies of the rapidly rotating secondary stars in 
CVs can also provide tests of stellar dynamo theories under 
extreme conditions. For example, questions regarding the 
impact of tidal forces on mag netic flux tube emergence (e.g. 
iHolzwarth fc Schiisslerl l2003h an d its effect on differential 
rota tion (e.g. Sgharlcmann 1982) are particularly pertinent 
fsee lWatson. Dhillon fc Shabaj|2006l for a discussion). 

Despite this, until recently there had been little di- 
rect observational evidence fo r magnetic activity in CVs. 
I Webb. Navlor fc Jeffries] (|2002l ) used TiO bands to infer the 
presence of spots on the secondary star in SS Cyg and es- 
timated a spot filling factor of 22 per cent. Unfortunately, 
this technique does not allow the surfaces of these stars to 
be imaged and hence the spot distributions could not be 
ascertained. 

Most recently. IWatson et al. I (|2006t ) used Roche tomog- 
raphy to map the starspot dist ribution on a CV sec ondary 
(AE Aqr) for the first time. In lWatson et all (120061 ) we es- 
timated that starspots covered approximately 18 per cent 
of the northern hemisphere of AE Aqr. The Roche tomo- 
gram of AE Aqr also showed that starspots were found at 
almost all latitudes, although there was a relative paucity 
of starspots at a latitude of ~40°. Furthermore, we found 
that, in common with Doppler images of single rapidly rotat- 
ing stars, AE Aqr also displayed a large high-latitude spot. 
In this work w e continue our series of papers on Roche to- 
mography fseelWatson fc Dhillonll200ll . IWatson et ai1l2003l . 
IWatson et aiT2 006 ) by mapping starspots on the long period 
(0.61-d) dwarf nova BV Cen for the first time. We also report 
on the serendipitous discovery of the binarity of HD220492. 



2 OBSERVATIONS AND REDUCTION 

Simultaneous spectroscopic and photometric observations 
were carried out over 3 nights on 2004 July 8-10. The spec- 
troscopic data were acquired using the 6.5-m Magellan Clay 
Telescope and the simultaneous photometry was carried out 
using the Carnegie Institution's Henrietta Swope 1.0-m Tele- 
scope. Both telescopes are situated at the Las Campanas 
Observatory in Chile. 



2.1 Spectroscopy 

The spectroscopic observations of BV Cen were carried 
out using the dual-beam Magellan Inamori Kyoc era Echelle 
spectrograph (MIKE - see iBernstein et al.|[20o3 ). The MIT 
Lincoln Labs CCD-20 chip with 2046 x 4096 pixels was 
used in the blue channel, and the SITe ST-002A chip, again 
with 2046 x 4096 pixels, was used in the red channel. The 



standard setup was used, allowing a wavelength coverage of 
3330A - 5070A in the blue arm and 4460A - 7270A in the red 
arm, with significant wavelength overlap between adjacent 
orders. With a slit width of 0.7 arcsec, a spectral resolution 
of around 38,100 (~7.8 km s" 1 ) and 31,500 (~9.5 km s" 1 ) 
was obtained in the blue and red channels, respectively. The 
chips were binned 2x2 resulting in a resolution element of ~ 
2.3 binned pixels in the red arm and 2.6 binned pixels in the 
blue arm for our chosen slit. The spectra were taken using 
400-s exposure times in order to minimise velocity smearing 
of the data due to the orbital motion of the secondary star. 
Comparison ThAr arc lamp exposures were taken every ~50 
minutes for the purposes of wavelength calibration. 

With this setup we obtained 63 usable spectra in each 
arm. Since the main goal of the Magellan run was to observe 
AE Aqr, we were restricted to 3-hour windows each night 
to observe BV Cen. Over our allocated 4 nights this would 
have allowed over 80 per cent of the orbit of BV Cen to be 
observed, but unfortunately we lost the final night due to 
bad weather. Other than that, the seeing was typically 0.6- 
0.7 arcsec on the first night and 0.9 arcsec on the next two 
nights, with occasional degradation to 1.5 arcsec. The peak 
signal-to-noise of the blue spectra ranged from 29 - 56 (typi- 
cally ~50) in the blue arm, and from 38 - 76 (typically ~65) 
in the red arm. Table Q] gives a journal of the observations. 

It should be noted that when using MIKE it is not possi- 
ble to change the slit orientation on the sky. In order to com- 
pensate for this, and reduce atmospheric dispersion across 
the slit, MIKE is mounted at a 30° angle to the Naysmith 
platform. This means that at zenith distances greater than 
~50° dispersion across the slit becomes significant. To avoid 
this, all observations of BV Cen were carried out at low air- 
masses and no exposure of BV Cen was carried out for an 
airmass above 1.39 (zenith distance > 44°). 

2.1.1 Data reduction 

The data were reduced using the MIKE redux IDL pipeline 
version 1.7. This automatically processes all calibration 
frames and then wavelength calibrates, sky subtracts, flux 
calibrates and optimally extracts the target frames. The fi- 
nal output consists of 1-d spectra for the blue and red arms. 
The typical rms scatter reported for the wavelength calibra- 
tion was around 0.002A. Since the redux package outputs 
vacuum wavelengths, whereas the linelists used in the Least 
Squares Deconvolution process (see Section [2} are measured 
in air, the wavelengths have been conv erted to air using the 
IAU standard given bv lMortonl (jl99l|). 

Unfortunately, we experienced some flux calibration 
problems around regions of strong emission lines (e.g. Ha 
and H/3), most likely due to poor tabulation of our flux stan- 
dard star HR 5501 around these lines. Since we mask out 
the strong emission lines during our study of the donor star 
absorption lines, this problem does not affect this work. We 
also found that we could not resolve (at the data reduction 
stage) a small jump in the flux level between the blue and 
red arms. We discuss the solution to this latter problem in 
Section [3] 

Since the secondary star contributes a variable amount 
to the total light of a CV, Roche tomography is forced to 
use relative line fluxes during the mapping process (i.e. it 
is not possible to employ the usual method of normalising 
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Table 1. Log of the spectroscopic observations of BV Cen, the relevant spectral- type template stars, a telluric-correction star and the 
spectrophotometric standards HR 5501 and HR 9087. The first column gives the object name, columns 2—4 list the UT start Date and 
the exposure start and end times, respectively, and columns 5-6 list the exposure times and number of spectra taken for each object. 
The final column indicates the type of science frame taken and, where applicable, the measured systemic velocities, 7, of the template 
stars computed from Gaussian fits to their Least Squares Deconvolved profiles (see Section [5}. 



Object UT Date UT Start UT End T exp (s) No. spectra Comments 



HR 5501 


2004 


Jul 


08 


22:40 


22:45 


1-7 


5 


Spectrophotometric standard 


BV Cen 


2004 


Jul 


08 


23:34 


02:31 


400 


22 


Target spectra 


HR 9087 


2004 


Jul 


09 


10:07 


10:13 


1-4 


6 


Spectrophotometric standard 


HR 5501 


2004 


Jul 


09 


22:40 


22:49 


1.5 


12 


Spectrophotometric standard 


BV Cen 


2004 


Jul 


09 


23:42 


02:21 


400 


20 


Target spectra 


Gl 863.3 


2004 


Jul 


10 


10:08 


10:14 


15-60 


4 


G5V template; 7 = +67.28 ± 0.10 km s" 1 


HD 221255 


2004 


Jul 


10 


10:25 


10:42 


20-100 


5 


G6V template; 7 = -11.55 ± 0.10 km s" 1 


HD 224287 


2004 


Jul 


10 


10:52 


11:04 


60-200 


5 


G8V template; 7 = +44.57 ± 0.10 km s~l 


HR 5501 


2004 


Jul 


10 


22:41 


22:47 


1.5-8 


8 


Spectrophotometric standard 


BV Cen 


2004 


Jul 


10 


23:35 


02:18 


400 


21 


Target spectra 


HD 217880 


2004 


Jul 


11 


10:20 


10:27 


45-80 


4 


G8IV template; 7 = -67.24 ± 0.10 km s" 1 


HD 220492 


2004 


Jul 


11 


10:32 


10:42 


45-120 


5 


G5IV template; binary 


BS 8998 


2004 


Jul 


11 


10:47 


10:52 


1.5-8 


6 


Telluric B9V star 



the data). In order to determine these relative fluxes, the 
slit losses need to be calibrated. The standard technique of 
using a comparison star on the slit is, however, not possible 
with MIKE due to the tightly-packed, cross-dispersed for- 
mat, and so simultaneous photometry is required (see Sec- 
tion [2]2]). We corrected for slit losses by dividing each BV 
Cen spectrum by the ratio of the flux in the spectrum (after 
integrating the spectrum over the appropriate photometric 
filter response function) to the corresponding photometric 
flux (see Section \2.2\ . We should note here that the slit-loss 
correction factors were only calculated using the red spec- 
tra, but applied to both arms of the spectroscopic data. A 
further correction was applied to the blue arm (to account 
for both a different slit-loss factor in the blue compared to 
the red and also the jump in flux level between the two arms 
mentioned earlier) at the Least Squares Deconvolution stage 
(see Section [3]). 

2.2 Photometry 

Simultaneous photometric observations were carried out us- 
ing a Harris V-band filter and the Site3 CCD chip with 2048 
x 3150 pixels. The CCD chip was windowed to a size of 
633x601 pixels (which covered the target and the brightest 
comparison stars, as well as the bias strip) resulting in a 4' x 
4' window on the sky. Windowing allowed the readout time 
to be reduced to 46-s for each 15-s exposure on BV Cen. 

2.2.1 Data reduction 

The photometry was reduced using standard techniques. 
Since the master bias-frame showed no ramp or large scale 
structures across the chip, the bias level of the frames was 
removed by subtracting the median of the pixels in the over- 
scan regions. Pixel-to-pixel sensitivity variations were then 
corrected by dividing the target frames through by a master 
twilight flat-field frame. Opti mal photometry was perform ed 
using the package photom (|Eaton, Draper fc Alien! I2OO2I ). 
There were two suitable comparison stars on each 
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Orbital Phase 

Figure 1. The V-band light curve of BV Cen. The data have 
been phased according to the ephemeris derived in Section [5] All 
3 nights data have been plotted together from phase 0-1. 



BV Cen target frame, identified from the ESO Guide 
Star Catalogue as GSC0866601471 (hereafter 1471) and 
GSC0866600859 (hereafter 0859). Two other bright stars 
in the field were deemed unsuitable since one was over- 
exposed, and the other had a nearby companion. Unfortu- 
nately, neither of the comparisons had reliable magnitude 
measurements available. We therefore used observations of 
the Landolt photometric standard SA105-437 to calibrate 
our instrumental magnitudes. This allowed us to determine 
the magnitude of 1471 and 0859 as m v = 13.43 ± 0.07 and 
m v — 13.54 ± 0.07, respectively. Finally, differential pho- 
tometry was performed using the bright nearby comparison 
1471. The light curve is shown in Fig. [1] 



2.3 Continuum fitting 

The continuum of the BV Cen spectra were fitted in ex- 
actly the same manner as that outlined in the Roche tomog- 
raphy study of AE Aqr (W atson et al.ll2006T l. Spline knots 
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were placed at locations in the spectra that were relatively 
line-free. The spline knot positions were chosen on an indi- 
vidual spectrum-by-spectrum basis (though the locations of 
the knots were generally the same for each spectrum) until a 
smooth and visually acceptable fit was obtained. Again, as in 
I Watson etail(|2006l ). we found that we systematically fit the 
continuum at too high a level, leading to continuum regions 
in the Least Squares Deconvolved (LSD - see Section[3]l pro- 
files lying below zero. The continuum fit was then shifted to 
lower levels until the continuum level in the LSD profiles lay 
at zero. The exact value of this shift varied from spectrum 
to spectrum, due to changes in the continuum shape as a 
result of variable accretion light, but the variation was gen- 
erally small. We found that this process did not change the 
actual LSD profile shape. 



3 LEAST SQUARES DECONVOLUTION 

Least Squares Deconvolution (LSD) effectively stacks the 
~1000's of photospheric absorption lines observable in a sin- 
gle echelle spectrum to produce a 'mean' profile of greatly 
increased signal-to-noise. The tech nique of LSD i s well doc- 
umented and was first applied by iDonati et al.l (|l997h and 
has since been used in many Doppler i maging studies (e.g . 



Barnes et alj 



Barnes et al.l 



19981; iLister et all 1 19991; iBarnes et all 120001; 



2001; Jcff ers et al.1 12002| ; IBarnes et all 12004; 



Marsden et al.l " 20050 as well as in the succe ssful mapping 



of starspots on AE Aqr (| Watson et al.|[2006h . For detailed 
informa tion on LSD, see these references as well as the re- 
view bv lCollier Cameronl l|200ll ). 

At present we use line lists gen erated by the Vi- 
enna Atomic Line Database (VALD - iKupka et al. 1 1 19991 ; 
iKupkaet ail fcoOO) for LSD. The spectral type of the sec- 
ondary star in BV Cen has been determined to li e within 
the range G5-G8 IV/V (|Vogt fe Brevsacherl 11980] ). In our 
analysis (see Section [5} later type spectral templates appear 
to fit the BV Cen spectra better, with our closest match 
being a G8IV template. In accordance with this, we down- 
loaded a line-list for a stellar atmosphere with T e / / = 5250 
K and log g — 3.55 (the closest approximation available in 
the database to a G8IV spectral type). We do not consider a 
slight error in the choice of line-list used as having a signifi- 
cant impact on the res ults of the LSD proc ess - see the dis- 
cussio n on this topic in lWatson et al.l (|2006l ) and also lBarned 
(1999). Since this line- list contained normalised line-depths, 
whilst Roche tomography uses continuum-subtracted data, 
each line-depth has been multiplied by the corresponding 
value in a continuum fit to the G8IV template spectrum. 

In order to circumvent the small flux 'jump' between 
the blue and red arms reported in Section [2.1. II we carried 
out LSD in each arm separately. For the red arm, LSD was 
carried out over a wavelength range of ~5050A to ~6860A, 
while for the blue arm the wavelength range was ~4030A to 
~5050A. Regions of strong emission lines such as Ha and 
H/3 were also masked out. Some 2790 lines were used in the 
deconvolution process across both arms, leading to a multi- 
plex gain in signal-to-noise of ~22 and yielding absorption 
profiles with a signal-to-noise of ~1300. Our version of LSD 
propagates the errors through the deconvolution process. 

Since the red-spectra were slit-loss corrected using si- 
multaneous photometry, the resulting LSD profiles from the 



red data are scaled correctly relative to one-another. We 
then rebinned the blue LSD profiles to the same velocity 
scale as the red profiles (since the blue arm had a superior 
resolution) using sinc-function interpolation. Slit-loss correc- 
tions to the blue LSD profiles were then made by scaling the 
blue LSD profiles to match the red profiles obtained at the 
same phase. This scaling was done using an optimal subtrac- 
tion algorithm. The blue profiles were scaled by a factor / 
and then subtracted from the red profile. The factor / that 
resulted in the minimum scatter in the residuals (compared 
to a smoothed version of the residuals) was then used to scale 
the blue LSD profiles. This method allowed the slit-losses in 
the blue to be calibrated despite the jump in flux between 
the red and blue arms. In addition, this method also allows 
further loss of light in the blue due to differential refraction 
to be corrected to first order. Deconvolving the blue and red 
data separately also allowed us to look for systematic prob- 
lems and/or differences in the profiles, of which none were 
apparent. We should note, however, that the LSD profiles 
exhibit a slight tilt, with the red continuum (more positive 
velocities) systematically higher than the continuum on the 
other side of the profile. This is a well-known artefact (e.g. 
lBarnedli"999h which will have little impact on the final image 
reconstruction. 

The individual LSD profiles are shown in Fig. [21 These 
clearly show the distinct emission bumps due to starspots 
moving from blue (-ve velocities) to red (+ve velocities) 
through the profile as BV Cen rotates. When the LSD pro- 
files are trailed (Fig. [3]) these features are still obvious in 
addition to the secondary stars orbital motion and varia- 
tions in the projected equatorial velocity, vsini. Fig. [3] also 
shows the trailed LSD profiles after removal of the orbital 
motion and subtraction of a theoretical profile at each phase. 
The theoretical profiles were calculated using our Roche to- 
mography code by adopting the binary parameters and limb 
darkening described in Section [6] and a featureless stellar 
disc. This process increases the contrast and enhances the 
starspot signatures (which now appear dark). 



4 ROCHE TOMOGRAPHY 

Roch e tomography is an alogous to Doppler Imaging 
(e.g. IVogt fc Penrodl I1983T I and has now been success- 
fully ap plied to the dono r star s of CVs on several oc- 
casions dRutten fc Dhillonl Il994l; iRutten fc Dhillonl Il996l; 
Watson et all 120031; ISchwope et al.1 |2004|; IWatson et all 



2006). Rather than repeat a detailed description of the 
methodology and axioms of Roche Tomography here, we 
refer the reader to the references above and the technical r e- 
views of Roche Tomograp hy by IWatson fc Dhillonl (12001 
iDhillon fc Watsonl l|200lli and IWatson fc Dhillonl i|2004 
The pertinent points with respect to this work are that we 
employ a moving uniform default map, where each element 
is set to the average value of the reconstructed map. Further- 
more, we d o not adopt a two-temperature or filling factor 
model (e.g. lCollier Cameron fc Unruh| [l994'). since these as- 
sume only two temperature components across the star while 
CV donors are expected to exhibit large temperature differ- 
ences due to the impact of irradiation. This means our Roche 
tomograms may be prone to the growth of bright pixels 
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Figure 2. The LSD profiles and maximum entropy fits (see Section[7]for further details) for BV Cen. The velocity of the secondary star 
has been removed using the parameters found in Section |6] and each profile is shifted vertically for clarity. The orbital phase for each 
exposure is indicated at the top right of each profile. 



(whic h two-temperature models suppress, ICollier Cameron] 
1 19921 ) and are quantitatively more difficult to analyse. 



5 EPHEMERIS 
5.1 Ephemeris 

Despite the brightness of BV Cen, and the fact that it is 
one of the longest period dwarf novae known, the system 
has been poorly studied. The most recent orb ital ephemeris 
for BV Cen was published bv lGillilandl 0982) and will have 
accumulated a large error over the last two decades. We have 
therefore determined a new ephemeris for BV Cen by cross- 
correlation with suitable template stars that were observed 
using the same instrumental setup. The cross-correlation 
was carried out over the spectral range 6400A - 6536A. This 
wavelength range not only covers several moderately strong 
lines and blends from the secondary that are useful for radial 
velocity measurements but include sev eral temperature and 
gravi ty sensitive lines for F-K stars (jStrassmci er fc Fekell 
ll99Ch . 



The BV Cen and template spectra were first rebinned 
onto the same velocity scale using sinc-function interpola- 
tion and then normalised by dividing through by a constant. 
The continuum was then subtracted by a 3rd order polyno- 
mial fit. This procedure is followed to ensure that the line 
strengths are preserved across the spectral region of inter- 
est. The template spectra were then broadened to account 
for the rotational velocity (y sin i) of the secondary star - 
the amount of broadening applied was determined using an 
optimal-subtraction technique. In this method the template 
spectra were first broadened by an arbitrary amount before 
being cross-correlated against the BV Cen spectra, allowing 
a first iteration on the radial velocity curve of BV Cen. The 
BV Cen data were then orbitally-corrected using the results 
of this radial velocity analysis and averaged. The template 
spectra were once more broadened by different amounts in 
0.1 km s _1 steps, multiplied by a constant, /, and then 
subtracted from the orbitally-corrected mean BV Cen spec- 
trum. The broadening that gave the minimum scatter in the 
residual optimally-subtracted spectrum was then applied to 
the template spectrum and the whole procedure repeated 
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Figure 3. Left-hand panel: A trail of the deconvolved profiles of BV Cen. The large gaps are due to non-continuous observations over 
three nights. The small gaps in the phase coverage are at times when arc spectra were taken for the purpose of wavelength calibration. 
Features due to starspots appear bright and several such features are clearly visible traversing the profiles from blue (more -ve velocities) 
to red (more +ve velocities) in the trailed spectra. Also evident is the orbital motion and variation in usini, which shows a maximum at 
phase 0.75 due to the varying aspect of the tidally distorted secondary star. Right-hand panel: The same as the left-hand panel, except 
the orbital motion has been removed using the binary parameters derived in Section [6] This allows the starspot tracks across the profiles 
to be more easily followed. In order to increase the contrast of this plot, we have subtracted a theoretical profile (sec text) from each 
LSD profile before trailing. Features due to starspots and irradiation appear dark in these plots. Note the narrow feature (indicated with 
an arrow and shown enlarged - see inset) that lies outside the blue edge of the profiles at phases 0.328-0.366, which seems to follow the 
motion of the large feature that runs through the profiles during this block of observations. This feature is discussed in more detail in 
Section 15] 



until the rotational broadening value no longer changed. 
This typically took 2 or 3 iterations. 

Through the above process a cross-correlation function 
(CCF) was calculated for each BV Cen spectrum. A ra- 
dial velocity curve can then be derived by fitting a sinusoid 
through the CCF peaks. Unfortunately, none of our spectral- 
type templates provided a good match to the BV Cen spec- 
tra, all requiring a multiplication factor / > 1, which would 
imply that the secondary star contributes more than 100 per 
cent of the system light. We also discovered that our G5IV 
template star HD220492 is actually a close binary (see Sec- 
tion Ek2J. 

Fortunately, as discussed in IWatson et al.l (|2006l ), the 
CCFs calculated are insensitive to the use of an ill-matching 
spectral template and, despite requiring multiplication fac- 
tors / > 1, the broadened template spectra still provided 
visually acceptable matches to the orbitally corrected BV 
Cen spectrum. The radial velocity curve in Fig. [4] was ob- 
tained after cross-correlation of the BV Cen spectrum with 
the G8IV star HD 217880 (our best guess spectral type for 
B V Cen) . This allowed us to obtain a new zero-point for the 
ephemeris of 



T =HJD 2453195.2859 ±0.0003 (1) 



with the orbit al period fixed at P = 0.611179 d (from 
lGillilandlll982r h This ephemeris was applied to the rest of 
the data presented in this paper. 

Cross-correlation with the G8IV template yielded a sec- 
ondary star radial velocity amplitude of K r = 137.3±0.3 km 
s _1 , a systemic velocity of 7 = -20.3 ± 0.2 km s _1 and a ro- 
tational broadening vsini = 95.3 km s _1 (see Section[6]for a 
discussion of the binary parameters) . These values were con- 
sistent within 2 km s _1 for all the spectral- type templates 
used. Since some of our spectral-type templates (including 
the G8IV template star HD 217880) did not have measured 
systemic velocities, we carried out least squares deconvo- 
lution (see Section |3} on each template star. The systemic 
velocity was then measured using a Gaussian fit to the sharp 
LSD profile, the results are listed in the final column of Ta- 
ble [1] We note that Gl 863 .3 has a published sys temic ve- 
locity (+66.9 ± 0.1 km s- MNordstrom et al.|[200j ') which is 
close to our measured value of +67.28 ± 0.10 km s _1 using 
this technique. 
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Figure 4. The measured secondary star radial-velocities with a 
sinusoidal fit from cross-correlation using a G8IV spectral-type 
template star. The data have been phase folded using the zero- 
point of Section [5] The error bars are typically smaller than the 
symbol size. 




Figure 5. The LSD profile of HD220492. This system is most 
likely a close binary system containing an early to mid G-type 
star and a later type companion. 



5.2 The binarity of HD220492 

While measuring the systemic velocities of the spectral-type 
templates, we discovered that the G5IV spectral-type tem- 
plate HD 220492 is a close binary system. The LSD profile 
of this system is shown in Fig.0and clearly shows two heav- 
ily broadened stellar line profiles. It is most likely that this 
is a close binary system containing an early to mid G-type 
star and later-type companion. We note that the binary na- 
ture of HD 220492 was missed by the Geneva-Copenhagen 
survey of the Solar neighbourhood, since no radial veloc- 
ity measurements of t his star appear to have been taken 
l|Nordstrom et alj2004l 1. HD 220492 is also listed as the opti- 
cal counterpart to th e X-ray source 1RXJS232426.0-450906 
i|Schwope et alj200Ch. and was i dentified as a suspected vari- 
able star by ISamus et al.l (|2004h . Given the large v sin i and 
hence rapid rotation of these stars, both components are 
likely to be magnetically active. Unfortunately, since only 
4 spectra over 7 minutes were taken, nothing can be said 
about the orbital period of this system. 



1 1.05 1.1 

secondary mass (M Q ) 

Figure 6. The entropy landscape for BV Cen, assuming an or- 
bital inclination of i = 53° and a systemic velocity of 7 = —22.3 
km s . Dark regions indicate component masses for which no fit 
to the target \ 2 could be found. The cross marks the position of 
maximum entropy, corresponding to component masses of Mi = 
1.18 M and M 2 = 1.05 M . 



6 SYSTEM PARAMETERS 

The system parameters (systemic velocity 7, orbital inclina- 
tion i, secondary star mass Mi and primary star mass Mi) of 
BV Cen h ave been determined using the same method as de- 
scribed in lWatson et alj {2006). In short, adopting the incor- 
rect parameters causes artefacts to appear in the Roche to- 
mograms leading to an increase in the structure (and hence 
a decrease in the entropy) of the final image. This can be 
visualised as an entropy landscape (see Fig. |5J), in which re- 
constructions to the same x 2 are carried out for different 
combinations of component masses and the entropy obtained 
for each set plotted on a grid of Mi versus M2. 

In order to search for the correct set of parameters we 
have constructed a series of entropy landscapes for different 
inclinations and systemic velocities. This is done in an iter- 
ative manner by first constructing entropy landscapes for a 
sequence of systemic velocities for a fixed inclination. The 
optimal systemic velocity found in this first iteration was 
then fixed and a series of entropy landscapes were then cal- 
culated out over a range of orbital inclinations. Once an op- 
timal orbital inclination was determined another sequence of 
entropy landscapes was carried out over a range of systemic 
velocities to ensure that the optimal systemi c velocity had 
not changed. As we have found previously l|Watson et al.l 
120031 ; IWatson et al.l (|2006T ) 'l . the systemic velocity obtained 
in this way is largely independent of the orbital inclination 
assumed in the reconstructions. 

For the reconstructions of BV Cen we have adopted a 
root-square limb-darkening law of the form 



= h [1 -a x (1-/0 - b x (1 - v^)] 



(2) 



where fi — cos 7 (7 is the angle between the line of sight and 
the emergent flux), 7o is the monochromatic specific inten- 
sity at the center of the stellar disc, and ax and bx are the 
limb-darkening coefficients at wavelength A. We calculated 
an effective central wavelength of 5067A for our spectro- 
scopic observations using 
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Figure 7. Points show the maximum entropy value obtained in 
each entropy landscape as a function of systemic velocity, assum- 
ing an orbital inclination of 53°. The solid curve is a 5th order 
polynomial which is only shown to emphasise the trend between 
entropy and systemic velocity. The highest data point corresponds 
to a systemic velocity of 7=— 22.3 km s _1 . 



(3) 



where di and at are the line-depths and error on the observed 
data at wavelength position Ai, respectively. This, therefore, 
takes into account the line depths and noise in the spectrum. 



Limb- darkening coefficients for a star of T e tf = 5250 K and 
log g = 3.5 (corresp onding to a G8 IV star - Pall et al]|2005l ) 
were obtained from Claret (1998) for the B- and V-bands. 
The values of the coefficients at these wavelengths were then 
linearly interpolated over to find a=0.5388 and fo=0.3088 for 
the effective central wavelength of our observations. These 
coefficients were used for all the reconstructions carried out 
for this paper. 

Fig.Qshows the peak entropy value obtained in entropy 
landscapes constructed assuming an orbital inclination of 
i=53° (this value was obtained after carrying out the iter- 
ative procedure described above) but varying the systemic 
velocity. This yielded an optimal value of 7 = -22.3 km s . 
We cannot assign a rigorous error estimate to the systemic 
velocity but we found that the image quality depended heav- 
ily on the assumed systemic velocity, and reconstructions 
were almost impossible for assumed systemic velocities that 
differed by more than ±2 km s _1 from the optimal value. In- 
deed, for this reason we sampled the systemi c velocity more 
finely than for AE Aqr l| Watson et alj|2006j ). incrementing 
in 0.1 km s _1 steps. It is unlikely, therefore, that the error 
on the systemic velocity exceeds ±1 km s _1 and is probably 
less than this. 

Our value of 7 = -22.3 km s , while close to the value 
of 7 = -20.3 ± 0.2 km s _1 derived from the radial veloc- 
ity analysis performed in Section [5] is in stark contrast to 
the two previo us published values known to us. These are - 
47±10 km s" 1 (|Vogt fc Brevsacherlll98Ch and -47±2 km s _1 
l|Gillilandi ri982). We are uncertain as to the cause of this dis- 
crepancy. It is unlikely that surface inhomogeneities, such as 
those caused by irradiation from the accretion regions (and 
which is known to cau se systematic errors in radial veloc- 

ity measurements - e.g. iDavev fc Smithlll992l ; IWatson et al.l landscapes vary as sin 3 i 



2003), could cause this discrepancy. First, our systemic ve- 
locity derived from radial velocity analysis only differs from 
that derived from the entropy landscape method by 2 km 
s~\ which implies that inhomogeneities on the surface of 
BV Cen have little impact . Second, one could argue that 
during the obse rvations of IVogt fc Brevsacheil (|l980T ) and 
iGillilandl l|l982h the impact of irradiation may have been 
much greater. In the case of HU Aqr, which is heavily ir- 
radiated, we found that the systemic velocity measured by 
standard radial velocity analyses c ould be shifted by a s much 
as -14 km s" 1 from the true value l|Watson et al.| [2003). Such 
a shift, if applic able to BV Cen, would bring t he results of 
IGillilandl (|l982l ) and IVogt fc Brevsacherl (|l980h into closer 
agreement with ours but only if a large portion of the lead- 
ing hemisphere (the side of the star as seen from <j) — 0.75) 
was irradiated during the ir observations. The lightcurve of 
IVogt fc Brevsacherl (1 19801 ) does indicate that the system was 
slightly brighter during their observations, but the system 
was certainly not in outburst and we would not expect a 
large amount of irradiation of the donor star. Furthermore, 
looking at our Roche tomograms in Section [7] irradiation 
of BV Cen appears to be fairly low, and certainly not large 
enough to impact a radial velocity curve analysis to the ex- 
tent needed to bring our values of the systemic velocity into 
agreement. Nor do we believe that our wavelength calibra- 
tion is incorrect, as we have measured the systemic velocity 
of Gl 863.3 using LSD and found it to be in agreement with 
iNordstrom et alj (120041 ') to within 0.4 km s -1 (see Section^. 
Therefore, we have no satisfactory explanation of the dis- 
cr epancy between our measu red s ystemic v e locity and those 
of IVogt fc Brevsacherl i| 19801 1 and IGillilandl (| 19821 ). 

Fig. [8] shows the maximum entropy values as a func- 
tion of the orbital inclination assuming the systemic veloc- 
ity of 7 = -22.3 km s" 1 derived earlier. This shows a clear 
trend with the best fitting inclination at i = 53° which we 
have adopted as the optimal value for BV Cen. This value 
is lower than (but agrees to within 2- a with) th e orbit al in- 
clinations of 62°± 5 ° and 61° ± 5° of IGillilandl (|l982h and 
IVogt fc Brevsacherl {l98Q), respectively. 

The entropy landscape for BV Cen carried out with i = 
53° and 7 = -22.3 km s _ is shown in Fig. [6] From this we 
derive optimal masses for the white dwarf and secondary star 
of Mi = 1.18 M and M 2 = 1.05 M , respectively. These 
values for the binary param eters seem to fall in b etwee n the 
previously published ones. IVogt fc Brevsacherl (|l980l ) also 
find relatively high masses for the binary components of 
Mi = 1.4 ± 0.2 and M 2 = 1.4 ± 0.2. Indeed, their de- 
rived white dwarf mass places it uncomfortably close to the 
Chandrasekhar limit while our new estimate lies more com- 
fortably within the allowed mas s limit for white dwarfs. On 
the other-hand, IGillilandl (Il982ft find lower masses of Mi = 
0.83 ±0.1 and M2 = 0.9 ± 0.1. These masses, however, lead 
to a high mass ratio (q = M 2 /Mi) of q = 1.08T0.18, which 
places BV Cen near the cr i tical mass ratio for mass tr ansfer 
instability (|Politandll99rj| ; iThoroughgood et al.ll2004T ). Our 
higher secondary star mass, and lower mass ratio (q — 0.89) 
moves BV Cen comfortably within the region where mass 
transfer is stable. We should note that, since the optimal 
inclination we find is quite low, a small error in the derived 
inclination leads to quite a large error in the masses. As 
expected, we find that the optimal masses in our entropy 



Roche tomography of CVs - IV. 



9 




Figure 8. Points show the maximum entropy value obtained in 
each entropy landscape for different inclinations, assuming a sys- 
temic velocity of 7 = -22.3 km s — 1 . The solid curve is a fifth or- 
der polynomial fit showing the general trend. The entropy value 
peaks at an orbital inclination of i=53°. Inclinations lower than 
49° result in an optimal white dwarf mass greater than the Chan- 
drasekhar limit and so we have not been extended the analysis to 
significantly lower inclinations. 



Figure 9. The Roche tomogram of BV Cen. Dark grey-scales 
indicate the presence of either star-spots or irradiated zones. The 
system has been plotted as the observer would see it at an orbital 
inclination of 53°, except for the central panel which shows the 
system as viewed from above the North Pole. The orbital phase 
(with respect to the ephemeris of Section [5]l is indicated above 
each panel. For clarity, the Roche tomograms are shown exclud- 
ing limb darkening. Individual spots referred to in the text are 
arrowed. 



Selection of the correct % 2 to aim for during the recon- 
structions is somewhat subjective. We select the x 2 which 
results in a close fit to the data (all reconstructions were 
carried out to x 2 =2.5, which indicates that our propagated 
errors are underestimated but this does not impact the final 
reconstructions), but not so close as to cause the Roche to- 
mograms to break up into noise. In order to determine how 
robust the derived parameters are, we have explored the 
effects that changing the aim \ 2 has on the entropy land- 
scapes. We find that fitting more closely to the data causes 
the maps to start to break up into noise, but we find that 
the derived inclination is the same (i — 53°), and the masses 
for each entropy landscape are typically within ±0.02 Mq, 
and we never see a difference greater than ±0.04 Mq for any 
assumed inclination. The same is true if we raise the aim \ 2 
(poorer fit to the data), except we lose sensitivity to the in- 
clination. Despite this, even fitting to substantially higher 
X 2 's where the maps become featureless we still find a pref- 
erence for inclinations between ~52 - 60°. We are therefore 
confident in the robustness of our derived binary parame- 
ters. 

Calculating the errors on our derived parameters is diffi- 
cult. One could perform a series of Monte Carlo simulations 
of the data and carry out entropy landscape reconstructions 
for each synthetic dataset. Unfortunately, this would require 
a large amount of time to complete. Instead we can deter- 
mine a rough guide to the errors from the scatter in the 
vsini and K r values obtained from conventional analysis of 
the data using the spectral-type templates described in Sec- 
tion [5] We found that the scatter in the measured K r and 
vsini values were, in both cases, ± 2 km s _1 . This leads to an 
error on the mass ratio, q, of ± 0.04 - which in turn trans- 
lates into an error on the component masses of ±0.03 Mq. 
Such an error agrees with the scatter we find in the optimal 
values returned from entropy landscapes constructed for dif- 
ferent x 2 ' s - This then gives optimal parameters of Mi = 1.18 



± 0.03 M and M 2 = 1.05 ± 0.03 Mq at an inclination of 
53°. 

Clearly, the errors in the masses are dominated by any 
error in the derived orbital inclination. Again, we can make 
a conservative estimate to the likely error in our orbital in- 
clination. A lower limit of i = 49° is set by the neccesity 
that the white dwarf mass be lower than the Chandrasekhar 
limit. If we assume that, therefore, our inclination is accu- 
rate to i=53±4°, we find Mi = 1.18 ±g;?§ M and M 2 = 
1.05 ±g;fl M . 



7 SURFACE MAPS 

Using the parameters derived in Section [6] we have con- 
structed a Roche tomogram of the secondary star in BV Cen 
(see Fig. [9}. The corresponding fit to the data are displayed 
in Fig. [2] A number of dark starspots are visible in the Roche 
tomogram. The most conspicuous of these is the high lat- 
itude spot (common in Doppler images of rapidly rotating 
stars) at a latitude of ~65° . Such high-latitude spots seen on 
rapidly rotating stars are commonly believed to be caused 
by strong Coriolis f orces acting to drive magn etic flux tubes 
towards the poles l|Schussler fc Solankilll992r ). This spot is 
als o similar to the large high-latitude spot seen on AE Aqr 
by IWatson et alJ (|2006h . not just in its latitude but also in 
the way that it appears towards the trailing hemisphere of 
the star. If such a shift in high latitude spots on CV donor 
stars is common then this hints that a mechanism is in op- 
eration that drives flux-tube emergence to this side of the 
star. We discuss this in more detail later. 

Also prominent is a dark feature near the L\ point. 
Although it is possible that this may be due to a large 
spot (with obvious consequences for theories that invoke 
starspots to quench mass transfer), we have to consider that 
the secondary stars in CVs are located close to an irradiating 
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source (the bright spot and/or white dwarf, for example). It 
is well known that the inner face of the secondary, and the L\ 
point in particular, is often irradiated. Indeed, early single- 
line Roche tomography studies of CVs concentrated on map- 
ping the so-called irradiation patterns on the inn er faces of 
the secondaries in CVs (e.g. IWatson et ah! [2003). While it 
is largely the contrast between the immaculate photosphere 
and lower spotted continuum contributions that cause spots 
to appear dark on our tomograms, irradiation causes the 
absorption lines to weaken due to ionisation and hence also 
appears dark. 

We believe that the feature near the L\ point in BV Cen 
is indeed due to irradiation. This interpr etation is supported 
by the modelling of B-V lightcurves bv lVogt fc Brevsacherl 
(1980) who found evidence for a slightly e nhanced t e mper - 
ature around the inner Lagrangian point. iGillilandl (|l982T l 
also found evidence for a slightly heated inner face. Close in- 
spection shows that the irradiation pattern is asymmetric, 
with the impact of irradiation strongest towards the lead- 
ing hemisphere. This is very similar to the irradiation pat- 
tern seen on the dwarf-no va IP Peg (see lWatson et al ] |2003l 
and lDavev fc SmitlJll992T) which has been explained by ir- 
radiation from the bright spot which is located on the cor- 
rect side of the star to cre ate this asymmetry. Certa inly the 
lightcurves of BV Cen (e.g. lVogt fe Brevsachedll980l) show a 
promi nent bright spot component. Also lVogt fc Brevsacherl 
( 1980) suggested that the B-V variations could be explained 
if a small region of the star nearest to the gas stream and 
bright spot had a temperature excess of AT 1 ~ 300 K. The lo- 
cation of this heated region described bv lVogt fc Brevsacherl 
(1980) matches what we observe in the Roche tomogram. 

Several other features in the Roche tomogram are also 
worth mentioning. One of these is a mid-latitude spot best 
seen at phase = 0.75 (marked as spot A on Fig[5J), which we 
have mentioned as it is one of the largest spots, other than 
the polar spot. Interestingly, there also seem to be two low- 
latitude spots near the L\ point (spots 'B' and 'C on Fig[5]), 
but located on the leading hemisphere. Although these spots 
are fairly heavily smeared in latitude, they may be located 
at a sufficiently low latitude to cross the mass-transfer noz- 
zle. It has been suggested that starspots are thought to be 
able to quench mass transfer as they pass the mass-losing 
nozzle on the don or, resulting in the low-states observed 
m m any CVs (see iLivio fc Pringlel 1199-4 iKing fc Cannizzol 
1998). The observations of low-latitude starspots near the 
L\ point lends some credence to these models - it seems ap- 
parent that sizeable starspots, which are often seen at high 
latitudes on rapidly-rotating stars, can also form at low lat- 
itudes near the L\ point in binaries. 

It is unlikely that the spots we have identified near the 
Li point are, instead, due to irradiation. This would im- 
ply that the irradiation is patchy in nature (rather than 
smoothly varying across the stellar surface) in order to form 
the dark spots we observe in the tomogram. This, in turn, 
would require that the irradiation is either beamed towards 
relatively small patches on the stellar surface, or would re- 
quire a complex accretion structure shielding the majority of 
the stellar surface and only allowing small spot-like regions 
to be irradiated. We can think of no reason why the irra- 
diation should be beamed in this manner, nor why such a 
complex accretion structure should exist. In addition, while 
the limb-darkening around the L\ region could be incorrect 




Figur e 10. Roche tomogra ms of the donor stars in AE Aqr (left 
- from IWatson et alj|2006l) and BV Cen (right) as viewed from 
above the pole (indicated with a cross). Both show a high latitude 
spot displaced towards the left-most (trailing) hemisphere (the 
orbital motion is towards the right). Furthermore, both appear 
to show a chain of spots descending from the polar spot towards 
the L\ point. 



by a significant degree due to the low effective-gravity, again 
this would be a smoothly varying artifact and would not lead 
to dark spots on the Roche tomogram. 

We also note that, when viewed from above the North 
pole of the star, there seems to be quite a distinct 'chain' 
of spots. This chain leads down from the polar regions to- 
wards the L\ point, showing a possible deflection towards 
the leading hemisphere with decreasing latitude. In Fig. [10] 
we pl ot an enlarged polar view of both BV Cen and AE Aqr 
(from IWatson et al.l 120061 ) which shows more clearly these 
features. Comparison with the Roche tomogram of AE Aqr 
shows a similar distribution of spots (we make a more de- 
tailed comparison in Section [9}. The fact that we see this 
same 'chain' of spots is suggestive of a mechanism that is 
forcing magnetic flux tubes to preferentially arise at these 
locations. Given that this is on the side of the star facing 
the white dwarf, this may be due to the impact of tidal 
forces which is thoug ht to be able to force spots t o arise at 
preferred longitudes (|Holzwarth fc Schu ssler 200 3j). In a re- 
cent D oppler Image of the pre-CV V471 Tau. iHussain et all 
(2006) also found the presence of high- latitude spots located 
on the side of the star facing the white dwarf. 

In addition to the features outlined above, a number 
of other small spots are visible. We have examined the in- 
tensity distribution of the pixels on our Roche tomogram 
in order to make a more quantitative analysis of the spot 
distribution. In our Roche tomography study of AE Aqr 
|Watson et al.| [2006). we looked for a bimodal distribution 
in pixel intensities and labelled the population of lower in- 
tensity pixels as spots, and higher intensity pixels as immac- 
ulate photosphere. Unfortunately, unlike AE Aqr, BV Cen 
does not show a clear bimodal distribution in pixel intensi- 
ties from which we can confidently distinguish between im- 
maculate and spotted photosphere (note that we neglected 
regions that are not visible and therefore have pixel intensi- 
ties set to the default value). Instead, the histogram of pixel 
intensities shows a broad peak, with long tails towards high 
and low pixel intensities. We have therefore defined a spot 
intensity by examining the polar spot feature. At the centre 
of the polar spot, the pixel intensities are very stable, and we 
have taken the intensity of the central regions of the polar 
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spot to represent the intensity of a 100% spotted region. Al- 
though there are lower intensity pixels present in the Roche 
tomogram (which contribute to the tail of low intensity pix- 
els in the maps), all of these were found to be confined to 
the irradiated zone near the L\ point. (We should note that 
the fact that the region around the L\ point has pixel in- 
tensities quite different from those of the spotted regions 
supports our interpretation that this feature is not a spot). 
We are therefore reasonably confident in our identification 
of the spotted regions on BV Cen. 

Determining the immaculate photosphere was some- 
what more difficult. There appear to be some bright regions 
in the map, especially near the polar spot. It is unlikely that 
these represent the immaculate photosphere; the growth of 
bright pixels in maps that are not 'threshol ded' is a known 
artifac t of Doppler imaging methods (e.g. lHatzes & Vogtl 
(|1992| ) - also see Section U), and it is likely that these fea- 
tures are due to this, contributing to the long tail of bright 
pixel intensities that we see. Instead, we have selected the 
intensity at the upper end of the broad peak of pixel inten- 
sities in the histogram to represent the immaculate photo- 
sphere. If we assume that the appearance and disappearance 
of spots simply alters the continuum level (and not the line- 
depths, which has a secondary influence on the LSD profiles) 
and a blackbody scaling, this gives a temperature contrast 
between the photosphere and spot of AT — 780K. Such a 
temperature difference seems reasonable. From simultane- 
ous modelling of lightcu rves and line-depth ratios of three 
active RS CVn systems, iFrasca et all (|2005h found a tem- 
perature difference between spotted an d immaculate photo- 
sphere of AT = 450-850 K. Similarly, iBiazzo et all (120061 ) 
find temperature differences of AT = 453-1012K for stars 
at different locations along the HR diagram, with stars of 
lower surface gravity having spots with a lower AT. 

Each pixel in our Roche tomogram was given a spot- 
filling factor between (immaculate) to 1 (totally spotted) 
depending on its intensity between our predefined immac- 
ulate and spotted photosphere intensities. After removal of 
the region near the L\ point (which is caused by irradiation 
and would cause us to overestimate the total spot cover- 
age), we find that some 25% of the visible hemisphere of BV 
Cen is spotted. This figure is relatively high compared to 
the spot coverages returned from Doppler images of isolated 
stars (typically 10 per cent) and is probably due to our less 
than ideal characterisation of spot filling-factors compared 
to imaging codes that invoke two-temperature models. How- 
ever, a visual compa rison of BV Cen with the map of AE Aqr 
(IWatson et al. 2006 - for which we estimated had 18% spot 
coverage) clearly shows that BV Cen has a higher coverage 
of spots on its surface, and our estimated spot coverage of 
25% for BV Cen is co n sisten t with this. Furthermore, a TiO 
study bv lWebb et ail (|2002T ) found a 22% spot coverage for 
the CV SS Cyg, comparable to what we find for BV Cen. 

Based on our spot characterisation scheme outlined 
above, we have determined the distribution of starspots as 
a function of latitude scaled either by the total surface area 
of the star, or by the surface area of the latitude strip over 
which the spot coverage was calculated (see Fig. Ill|l . This 
shows that the highest spot-filling factor is achieved at polar 
latitudes, primarily due to the large high-latitude spot cen- 
tred at +65°. We also find that the spot coverage reaches a 
minimum at intermediate latitudes around ~45° . This is also 



very similar to what was found for AE Aqr l|Watson et al.l 
2006), though the drop in spot coverage at this latitude is 
far more evident in BV Cen than for AE Aqr. 

Fig. [TT] also shows clear evidence for a bimodal distri- 
bution of starspots with latitude. Our analysis indicates a 
lower-latitude site of spot emergence around 25° ; indeed the 
majority of spots seem to form at lower latitu des in BV Cen. 
In ou r Roche tomography study of AE Aqr l|Watson et al.l 
2006), we found some evidence for spots at lower latitudes, 
although this was somewhat uncertain given the low signal- 
to-noise of those observations. While the BV Cen data has 
large phase gaps, this time we can be confident that we 
are seeing spots at lower latitudes. Many of the features 
described here, such as the off-pole high latitude spot, the 
bimodal spot distribution and the relative paucity of spots 
at intermediate latitudes are also seen in Dop pler Images of 
a num ber of other stars (see the discussion in I Watson et alj 
120061 and references therein). What is most startling, how- 
ever, is the great similarity betw een the maps of AE Aqr 
presented in I Watson et ai1 {2006), and those of BV Cen pre- 
sented in this work. 

Finally, we have checked whether the features seen in 
the Roche tomogram are real, or artefacts due to noise. We 
carried out two further Roche tomography reconstructions, 
the first only using the odd numbered spectra, and the other 
using the even numbered spectra. These are both shown in 
Fig. [12] and, despite the reduction in phase coverage, both 
maps show the same features. We are therefore confident 
that the features in the Roche tomogram are not due to 
noise. 



8 A SLINGSHOT PROMINENCE? 

In addition to the many starspot signatures that are visible 
in BV Cen's trailed spectra (see Fig|3}, a curious narrow fea- 
ture is also evident between phases 0.328 and 0.366 on the 
blue edge of the profile. This feature appears as a narrow 
continuation of the main track through the profiles during 
this block of observations, but lies outside the stellar absorp- 
tion profile and therefore cannot lie on the stellar surface. 
Closer inspection of the individual profiles (see Fig. I13p re- 
veals a narrow, weak emission feature in the continuum of 
5 LSD profiles. In order to confirm the reality of this fea- 
ture, we have visually inspected the profiles that were de- 
convolved from the blue and red spectra separately to see 
whether they are present in both sets of profiles. Indeed this 
feature is present (albeit very weakly) in both sets of profiles 
which makes it unlikely to be due to noise or a systematic 
effect arising during the LSD process. This feature is also 
not due to contamination from lunar or solar light during 
the obser vations as this would appear as an absorption fea- 
ture (e.g. iMarsden et al.ll2005h . Furthermore, moon-rise did 
not occur until at least 2 hours (extending to 4 hours on the 
final night) after observations of BV Cen were concluded. 
The fact that the feature is at zero velocity with respect to 
BV Cen's systemic velocity of -22 km s _1 rules out a 'ter- 
restrial' origin, which would be centred at km s~ . We are 
therefore confident that this feature is real. 

Since the emission feature lies outside the stellar line 
profile, and hence lies off the stellar limb, it can only be 
attributed to circumstellar material. Solar prominences ap- 
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Figure 11. Plots showing the spot coverage on BV Cen as a function of latitude. Left: spot coverage as a function of latitude expressed 
in terms of the surface area at that latitude. Right: spot coverage as a function of latitude normalised by the total surface area of the 
northern hemisphere. Since the grid in our Roche tomograms is not aligned along strips of constant latitude, some interpolation between 
grid elements is required in order to produce these plots. This results in their slightly noisy nature. 




Figure 12. As for Fig. [9] but reconstructions for odd-numbered spectra (left) and even-numbered spectra (right). 



pear as bright emission loops when they are viewed off the 
solar limb, and it is probable that we are also seeing a 
prominence structure on BV Cen. Indeed, large prominences 
have been re ported on other rapidly rota ting stars such as 
AB Dor (e.g. ICameron fc Robinsonlll989l ) and Speedy Mic 
ijDunstone et al.ll2006al ) and are often observed as transient 
absorption features passing through the Doppler broadened 
Ha stellar line. Recently, howev e r, anal ysis of VLT data of 
Speedy Mic by iDunstone et alj |2006bl ) has revealed rota- 
tionally modulated emission outside of the stellar Ha line 
due to loops of emission seen off of the stellar disc, but which 
can also be associated with prominences seen to transit the 
stellar disc at other times. In addition, peculiar low-velocity 
emission features seen in SS Cyg and IP Peg during out- 
burst have also been interpreted as 's l ingsho t prominences' 
jSteeghs et al.lll99fj) . |Gaensicke et all (fl998) discovered the 
presence of highly-ionized low velocity-dispersion material 
located between the L\ point and the centre of mass in AM 
Her which they attributed to a slingshot prominence. Sim- 
ilarly, triple-peaked Ha lines following the m otion of the 
donor star in AM Her has been reported b y IKafka et al.l 
|2005l ) and lKafka, Honevcutt fc Howelll (|2006h - The authors 



interpreted these as long-lived prominences on the donor 
star and noted that one component was consistent with stel- 
lar activity lying vertically above the L\ point. 

The emission feature we see in BV Cen appears sta- 
tionary at ~0 km s _1 within the binary frame. This is in 
keeping with the statio nary slingshot prom inences seen in 
SS Cyg and IP Peg by ISteeghs et all il996T). and the low- 
veloc ity emission observed in other CVs ( Marsh fc Home! 
1990). Generally, we would expect the emission from promi- 
nences observed off the stellar limb to be weak and unde- 
tectable for CVs given their faintness. However, it is possible 
that the prominences are illuminated by the accretion light, 
causing prominences forming between the donor star and the 
white dwarf to become visible when otherwise they would 
be undetectable. Certainly, the low velocity of the emission 
suggests a position close to the centre-of-mass of the binary 
at a point between the donor star and white dwarf where 
such illumination is most likely. 

Given that prominences are normally only seen in lines 
that form above the photosphere (e.g. the Hydrogen Balmer 
lines), observing them in the LSD profiles which are obtained 
from photospheric lines is unexpected under normal condi- 
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tions. We believe that the emission seen in BV Cen's photo- 
spheric lines is due to excitation of these species within low 
density gas in the prominence due to the impact of irradia- 
tion. Thus, we suspect that not only does irradiation cause 
the prominence to become highly visible, but also causes it 
to be observable in some photospheric lines in which promi- 
nences in a normal, unirradiated environment would not nor- 
mally emit. We have considered that the emission could be 
due to a wind launched from the accretion regions. Such a 
wind, however, would exhibit a radial velocity modulation 
due to the orbital motion of the primary star which we do 
not observe. Furthermore, presumably the velocities of ma- 
terial in such a wind would produce a far broader emission 
line than observed in BV Cen. For these reasons, and the 
fact that it kinematically matches previous observations of 
slingshot prominences in CVs, we prefer the interpretation 
that this feature is the result of an irradiated prominence. 

8.1 Limits on the prominence size and height 

We find that the emission feature in BV Cen is very narrow, 
with a velocity width ( AV) of ~ 10 km s -1 . Using this width 
we can pl ace an upper l i mit o n the emission source size, I. 
Following ISteeghs et al.l (|l996h , we assume the prominence 
is co-rotating with the secondary star, so we can write, 



K x +K 2 a yj 

where Kj and K% are the radial velocity amplitudes of the 
white dwarf and donor star, respectively, and a is the orbital 
separation. For the parameters derived for BV Cen in Sec- 
tion[6l this places an upper limit of 75,000 km. Naturally, the 
velocity dispersion is increased by instrumental resolution, 
thermal broadening, saturation broadening and turbulence 
within the prominence itself. For a 10,000 K prominence 
the thermal Doppler veloci ty of Hydrogen is 12.9 km s^ 1 
(e.g. iDunstone etalll2006bh . Since the lines included in our 
LSD are mainly heavier elements such as Fe and Ca, the 
thermal broadening will be far smaller at ~2 km s . Solar 
pro minences exhibit turbule nt motions of several km s , 
and IDunstone et alj l|2006bl ) estimate a turbulent velocity 
of 5 km s _1 for the prominences observed on Speedy Mic. 
However, the dominant broadening mechanism is the ~9.5 
km s _1 instrumental resolution of our observations. Consid- 
ering these limitations, our estimated maximum source size 
of 75,000 km should be viewed with caution as it is possible 
that the source is actually unresolved. 

We have further analysed the LSD profiles to check 
whether or not we can see the zero-velocity prominence as 
it tracks across the stellar disc. We were unable to posi- 
tively identify any feature with zero- velocity between phases 
0.374-0.522 during the first night's data. This may be due 
to the large emission bump feature that traverses the profile 
due to irradiation and/or starspots (see Section [7]) which 
probably masks the prominence itself. Analysis of the sec- 
ond night's data, however, does reveal a narrow emission 
bump at zero-velocity between orbital phases 1.974-2.058. 
This is quite conspicuous (see Fig. 1 13 j> as the emission fea- 
ture moves across the profile in the opposite direction to the 
starspot features, and Roche Tomography is clearly unable 
to fit it. This agrees with the picture of a prominence holding 
material near the centre-of-mass of the binary, and that at 



orbital phase ~0 we are effectively looking over the top (or 
possibly under the bottom) of the star at the prominence on 
the other side. This, therefore, means that the prominence 
structure must be raised above or below the orbital plane 
in order for it not to be eclipsed by the donor star at this 
phase. 

Given the lack of eclipse and the parameters found for 
BV Cen in Section[6] combined with the assumption that the 
prominence is located above the centre-of-mass of the binary 
(the point of zero- velocity) , we find that the prominence 
must lie at least 160,000 km above the orbital plane. If we 
assume that the prominence lies below the centre-of-mass 
(i.e. we are looking 'underneath' the star around phase 0) 
then it lies out of the orbital plane by at least 2,400,000 
km. We feel that this latter case is highly unlikely. First, a 
prominence below the orbital plane is far more likely to be 
eclipsed by the accretion regions around orbital phase ~0.35, 
when it is clearly seen in the data. Second, if we assume that 
the prominence is only visible due to illumination from the 
accretion regions then a prominence below the orbital plane 
is located too far away from the irradiating source. 

8.2 Prominence evolution and structure 

Three nights of consecutive (albeit interrupted) observations 
allows a limited discussion of the evolution of the slingshot 
prominence observed on BV Cen. The prominence is cer- 
tainly seen at the start of the first night (<f> — 0.328-0.366) 
before it transits the stellar disc, where it then becomes in- 
visible. We have assumed that, rather than the prominence 
disappearing at this point, its signature is lost in the com- 
plex structure present in the stellar line profile. We are able 
to pick up an emission feature with the same position in 
velocity space again at the start of the second night (<f> = 
1.974-2.038, see Fig.flSJ). 

Curiously, this prominence feature then disappears after 
orbital phase 2.038. Since, unlike the first night, the promi- 
nence on the second night appears quite clearly in the middle 
of the stellar line profile it is, on this occasion, difficult to 
explain how the feature could suddenly be lost within the 
stellar line. This suggests that the disappearance may be due 
to rapid evolution of the prominence itself. Indeed, it does 
seem as though the emission feature weakens during the sec- 
ond night before its apparent disappearance (see Fig. 113ft . 
Though the interpretation of a prominenc e rapidly evolv- 
ing on timescales of hours is speculative, IDunstone et al.l 
(2006b) also find evidence of individual prominences evolv- 
ing on timescales of ~9 hours on Speedy Mic. Certainly, 
we can find no evidence for the prominence feature on the 
3rd night. This is despite the fact that the emission should 
be well separated from the stellar line profile after orbital 
phase 3.645. This supports the idea that the prominence ma- 
terial is not long-lived, only lasting a couple of days before 
either draining back to the stellar surface or being ejected. 
IDunstone et al.l (|2006al ) found that, while some prominences 
on Speedy Mic were still visible after 5 nights, others formed 
or disappeared over the course of one night. 

Finally, unlike the prominences seen on Speedy Mic and 
AB Dor which lie typically between 2 to 9 stellar radii from 
the stellar rotation axis (with a concentration at the co- 
rotation radius), the feature observed in BV Cen is far closer 
to the stellar surface. If the prominence is located above the 
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centre-of-mass of the binary, then this places it just 1.5 R* 
(assuming the Roche-lobe volume radius for BV Cen) from 
the rotation axis of the secondary star. Although this is at 
odds with most Ha observations which show prominences 
at or beyond the co-rotation radius, clouds substantially 
closer to the stellar surface hav e been reported in HK Aqr 
l|Bvrne. Eibe fc Rollestonlll996l reported prominences 0.34 - 
3.2 R* above the s tellar surfac e) and RE J1816+542 (heights 
as low as 0.88 R*. lEibdli~998h . 

The height of the prominence feature seen in BV Cen 
also agrees with X-ray observations of rapidly-rotating iso- 
lated stars, even when Ha observations of prominences of 
the stars in question reveal high cloud heights . For instance, 
in Ch andra X-ray observations of AB Dor, iHussain et al.l 
l|2005l ) found that a significant fraction of the emission arose 
from compact regions near the stellar surface with heights of 
less than 0.3 R», and that the emitting corona does not ex- 
tend more than 0.75 R, above the stellar surface. Obviously, 
the formation mechanism of these prominences and active 
regions necessitates more work. In particular, the apparent 
preference for 'slingshot prominences' observed in CVs to 
form at low-velocity sites requires satisfactory explanation. 



9 DISCUSSION 

We have shown that the donor star in BV Cen is highly spot- 
ted, the second CV donor for which this has been shown 
to be the case. Fig. [10] highlights the striking similarities 
in spot distribution between BV Cen and AE Aqr, despite 
their quite different binary parameters and spectral types. 
Although BV Cen has a higher spot coverage and more low- 
latitude spots than AE Aqr, both show a high latitude spot 
displaced towards the trailing hemisphere. It would be inter- 
esting if all donor stars in CVs showed high latitude spots 
displaced in the same direction. Such a deflection could be 
the result of the orbital motion of the binary. Since the ro- 
tation axis lies outside the donor stars in these binaries, one 
may expect the play of Coriolis (and/or centrifugal) forces 
on magnetic flux tube emergence to be different from sin- 
gle stars. Further observations, however, are required before 
we can definitely say that such a deflection of high-latitude 
spots takes place in CV donors. 

Another striking similarity between the two tomograms 
presented in Fig. QJj] is the apparent chain of spots extend- 
ing down from the polar regions to the L\ point. This is 
more evident in the tomogram of AE Aqr, but probably only 
because there are fewer low latitude features which allows 
these spots to stand-out more on AE Aqr. We believe that 
this chain of spots is probably due to the impact of tidal 
forces due to the close proximity of a compact companion, 
and such a 'sub-white dwarf c oncentration of spo ts is also 
seen on the pre- CV V471 Tau (|Hussain et al.ll200rj ). A con- 
centration of starspots on the inner face of the donor stars in 
CVs may also have consequences for the accretion dynamics 
of these objects. It has long been thought that starspots may 
be able to quench mass-transfer from the donor as they pass 
across th e L\ point, leading to the low-states seen in many 
CVs (e.g. lLivio fc Pringlei[l994l . lKing fc Cannizzoll 19981 ). In- 
deed, in their study of the ma s s-tran sfer history of AM Her, 
iHessman. Gansicke fc Mattel ([20001 ) concluded that such a 
model would require an unusually high spot-coverage near 
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Figure 13. Plots of the LSD profiles that show the narrow emis- 
sion feature at a velocity of ~0 km s _1 on nights 1 & 2 as dis- 
cussed in the text. The Roche tomography fits to the data are 
also plotted and the orbital phase is indicated on the left. Unlike 
Fig. [2] the orbital motion has not been removed and a vertical line 
centred on a radial velocity of km s — 1 has been plotted to show 
where the emission feature appears. The slingshot prominence 
feature appears outside of the stellar lines on the first night's 
data (top 6 profiles). During the second night (bottom profiles) 
it appears within the profiles but moves in the opposite direction 
relative to the starspot features and the Roche tomography code 
is therefore unable to fit this feature. 
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the Li point, or otherwise some mechanism that drives spots 
towards the L\ point. Certainly, both AE Aqr and BV Cen 
do seem to show increased spot coverages towards the L\ 
point in support of their conclusions. 

The fact that we see more active regions near the L\ 
point may also explain why we appear to see a preference 
for 's lingshot prominences' to form above the donors inner 
face (Steeghs ct al. 1996; Gacnsicke et al. I I 19981 ; iKafka" et al.l 
120051 ; IKafka et al.ll2006l - and again in this work) . Certainly, 
the fact that prominence material at these locations will be 
illuminated by the accretion regions also means that ob- 
servations will be biased towards detecting prominences in 
the region between the white dwarf and donor star in CVs. 
Furthermore, if surface magnetic fields are strong enough in 
the neighbourhood of the L\ point this may also cause frag- 
mentation of the mass flow, resulting in the inhomogeneous 
or 'blobby' accretion se en in some CVs (e.g. lMeintiesll2004 
iMeinties fe Juruall2006h . 



10 CONCLUSIONS 

BV Cen is the second CV donor star for which starspots 
have unambiguously been imaged on. Again, as with AE 
Aqr, we find a high (25 per cent) spot coverage, and the high 
activity-level is further confirmed by the detection of a sling- 
shot prominen ce. Compariso n with the Roche tomograms of 
AE Aqr from I Watson et all |2006) show many similarities 
between the two systems despite the quite different funda- 
mental parameters of these two binaries. These spot distri- 
butions hint at the impact of tidal and/or Coriolis forces on 
the emergence of magnetic flux tubes in these binaries, and 
suggest that the inner faces of CV donor stars are unusually 
heavily spotted. As always, further observations are recom- 
mended in order to confirm that such spot distributions are 
widely seen on CV donors. If a fixed spot distribution were 
to be found, this would also suggest that differential rotation 
is suppressed (or at lea st weak) in CV donors, as suggested 
by IScharlemannl (j 19821 ) . 
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